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I.  Introduction 


The  Defense  Nuclear  Agency,  through  the  Army  Corps  of  Engineers,  is  in  the  process 
of  designing  a  facility  called  the  Large  Blast/Thermal  Simulator  (LB/TS).  The  US  Army 
Ballistic  Research  Laboratory  (BRL)  has  been  assigned  the  task  of  developing  the  blast 
and  thermal  simulation  technology  necessary  for  the  design  of  this  facility.  The  facility  is 
intended  to  subject  military  vehicles  such  as  tanks  and  helicopters  to  the  blast  and  thermal 
loading  associated  with  a  nuclear  burst.  An  illustration  of  the  proposed  facility  can  be  seen 
in  Figure  1. 


The  LB/TS  is  essentially  a  large  shock  tube  in  which  compressed  gas  is  quickly  released 
from  a  reservoir  called  the  driver  to  form  a  shock  wave  which  travels  down  the  expansion  sec¬ 
tion  of  the  shock  tube.  By  adjusting  the  volume  of  the  driver  and  the  initial  driver  pressure, 
a  desired  combination  of  shock  overpressure  and  simulated  weapon  yield  can  be  obtained. 
Thermal  radiation  is  simulated  through  combustion  of  an  aluminum/oxygen  mixture  near 
the  target. 

In  shock  tube  operation,  when  a  shock  wave  meets  the  open  end  of  the  expansion  section 
a  rarefaction  wave  is  generated.  This  rarefaction  wave  travels  back  into  the  expansion  section 
in  the  opposite  direction  of  the  initial  shock  wave  and  degrades  the  flow  characteristics  of 
the  shock  wave,  thus  ruining  the  nuclear  blast  simulation. 

In  a  small  shock  tube  the  obvious  solution  to  this  problem  is  to  make  the  expansion 
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section  long  enough  that  the  rarefaction  wave  does  not  arrive  at  the  test  section  until  after 
the  period  of  interest  has  past.  This,  however,  is  not  a  feasible  alternative  for  a  facility  as 
large  as  the  LB/TS.  As  currently  conceived,  the  expansion  section  of  the  LB/TS  has  a  cross 
sectional  area  of  163  m2  and  a  length  of  160  m.  Computational  studies  have  shown  that 
a  design  without  an  RWE  would  require  a  much  longer  expansion  section.  For  a  13.4  kPa , 
600  KT  simulation  an  expansion  section  length  greater  than  1500  m  would  be  required. 
Obviously  the  cost  of  building  such  a  facility  would  be  prohibitive. 

The  alternative  solution  to  this  problem  is  to  use  a  device  called  an  active  Rarefaction 
Wave  Eliminator  (RWE)  which  would  be  attached  to  the  end  of  the  expansion  section.  An 
RWE  is  essentially  a  converging  nozzle  whose  exit  plane  cross-sectional  area  can  be  changed 
with  time.  As  the  flow  passes  through  the  RWE  it  is  accelerated  by  the  area  reduction  so 
that  the  static  pressure  of  the  flow  is  equal  to  that  of  the  atmosphere  or  the  Mach  number 
at  the  RWE  exit  is  equal  to  1.  If  one  of  these  two  conditions  exists  then  no  rarefaction  wave 
is  generated.  The  RWE  exit  plane  area  is  continuously  changing  to  satisfy  these  conditions 
for  the  varying  flow. 


II.  The  1:57  Scale  Shock  Tube 

To  verify  the  theory  of  RWE  operation  developed  at  the  University  of  Toronto  and 

1,  2,  3 

BRL  ,  a  small  scale  active  RWE  was  installed  on  a  shock  tube  at  BRL.  This  shock 
tube,  shown  in  Figure  2,  has  an  expansion  section  diameter  of  25.40  cm,  a  single  cylindrical 
driver  whose  diameter  is  10.16  cm  and  a  throat  diameter  of  6.40  cm.  When  compared  to  the 
proposed  full  scale  LB/TS  this  shock  tube  scales  linearly  to  1:57  of  its  size.  All  diameters 
noted  here  and  hereinafter  are  internal  diameters  unless  stated  otherwise. 

The  shock  is  initiated  by  the  natural  burst  of  a  diaphragm  located  in  the  throat.  The  de¬ 
sired  driver  pressure  determines  the  appropriate  thickness  of  diaphragm.  During  the  testing, 
three  different  lengths  of  driver  section  were  used.  These  lengths  were  30.48  cm,  55.88  cm 
and  106.68  cm  resulting  in  driver  volumes  of  2840  cm3,  4988  cm3  and  9018  cm3  respectively. 

Figure  2  shows  that  the  expansion  section  length  was  set  at  either  856.61  cm  or 
1719.58  cm  during  the  tests.  The  longer  expansion  section  was  used  to  generate  baseline 
data  which  had  no  rarefaction  wave  during  the  positive  phase  of  the  flow.  This  baseline  data 
served  two  purposes.  Primarily,  the  data  from  the  long  expansion  section  tests  was  used 
to  generate  the  driving  functions  for  the  RWE.  Also,  the  pressure  histories  obtained  from 
the  long  expansion  section  tests  were  the  standard  by  which  the  RWE  tests  were  judged. 
The  shorter  expansion  section  was  used  during  tests  with  and  without  the  RWE.  The  short 
expansion  section  tests  with  no  RWE  were  performed  to  show  what  effect  the  rarefaction 
wave  had  on  the  pressure  histories.  This  same  expansion  section  length  was  used  in  the  tests 
of  the  RWE. 

A  total  of  nine  channels  of  data  is  recorded  for  each  test  of  this  shock  tube.  A  static 
pressure  transducer  is  used  to  measure  the  pressure  in  the  driver.  A  stagnation  probe  is 
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1.  ALL  DIAMETERS  INSIDE 

2.  DIMENSIONS  IN  CENTIMETERS 

3.  NOT  TO  SCALE 

4.  RWE  NOT  SHOWN 


Figure  2.  The  25.40  cm  Diameter  Shock  Tube 

located  in  the  expansion  section  52.1  cm  downstream  from  the  diaphragm.  This  gage  is  used 
to  trigger  the  recording  sequence  and  the  EWE.  A  pair  of  static  and  stagnation  probes  is 
mounted  in  the  expansion  section  at  three  locations,  resulting  in  six  additional  gages.  The 
first  pair  of  gages  is  located  7.519  m  from  the  diaphragm.  The  second  and  third  pairs  are  at 
7.900  m  and  8.281  m,  respectively.  Finally,  a  rotary  potentiometer  was  used  to  measure  the 
angle  of  attack  of  the  RWE  louvers.  The  RWE  resembles  a  shutter  valve  and  the  angle  of 
its  louvers  determines  the  amount  of  area  available  for  the  gas  to  flow  through  the  device. 
Thus  the  rotary  potentiometer  makes  it  possible  to  record  the  exit  plane  area  time  histories. 

The  location  of  the  three  static  and  stagnation  pairs  is  significant.  The  test  station 
located  8.281  m  from  the  diaphragm  is  the  station  closest  to  the  eventual  placement  of  the 
RWE.  Therefore  the  data  collected  at  this  station  is  used  to  generate  the  driving  functions 
for  the  RWE.  The  station  at  7.519  m  is  important  because  the  distance  from  this  station  to 
the  RWE,  when  scaled,  is  equivalent  to  the  distance  from  the  downstream  end  of  the  test 
section  to  the  RWE  in  the  full  scale  LB/TS.  The  data  collected  at  this  station  ultimately 
determines  whether  or  not  the  RWE  functioned  successfully.  The  remaining  station  is  to  be 
used  as  a  backup  in  the  event  that  gages  at  one  of  the  other  stations  fails  during  a  test. 
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III.  The  1:57  Scale  Active  Rarefaction  Wave 

Eliminator 

The  active  rarefaction  wave  eliminator  which  was  tested  is  illustrated  in  Figures  3  and 
4.  In  section  I  the  RWE  was  described  as  a  simple  converging  nozzle  placed  at  the  open 
end  of  a  shock  tube.  Inspection  of  these  drawings  shows  that  in  reality  the  RWE  is  much 
more  complex.  While  it  is  possible  to  computationally  model  a  conical  nozzle  whose  shape 
changes  with  time,  designing  a  machine  to  do  the  same  is  quite  difficult.  The  end  result 
is  the  device  shown  in  the  figures.  This  RWE  uses  rotating  louvers  to  block  a  controlled 
portion  of  the  exit  area  thus  creating  the  same  effect  as  a  converging  nozzle. 


DIMENSIONS  IN  CENTIMETERS 


Figure  3.  The  1:57  Scale  RWE  (front  view) 


1.  Louver  Overlap 

By  examining  the  illustrations,  one  can  see  that  as  a  louver  rotates,  the  amount  of  its 
area  which  is  presented  to  the  shock  tube  is  proportional  to  the  sine  of  the  angle  between 
the  louver  and  the  axis  of  the  shock  tube.  Consider  for  a  moment  a  single,  thin,  flat  plate 
in  space.  When  the  plate  is  nearly  horizontal,  a  modest  angular  displacement  will  result 
in  a  sizable  change  in  area  presented  to  a  flow  field  traveling  horizontally.  However,  when 
the  plate  is  nearly  vertical  a  great  change  in  angle  is  required  to  achieve  the  same  change 
in  presented  area.  From  this  it  can  be  seen  that  as  the  attack  angle  of  the  plate  increases, 
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DIMENSIONS  IN  CENTIMETERS 

Figure  4.  The  1:57  Scale  RWE  (side  section) 

its  angular  velocity  must  increase  rapidly  for  a  constant  rate  of  area  change.  This  high 
acceleration  rate  places  a  large  torque  requirement  on  the  motor  and  the  rotating  louvers. 

To  overcome  this  problem  the  louvers  were  designed  in  such  a  way  that  they  overlap. 
Adjacent  louvers  are  offset  by  an  angle  of  5  degrees  so  that  they  will  not  interfere  with  one 
another.  As  a  result,  the  overlap  first  occurs  at  a  louver  angle  of  about  61  degrees.  When  the 
overlap  condition  occurs,  the  area  of  the  shock  tube  is  completely  blocked  by  the  projected 
area  of  the  louvers.  This  design  greatly  reduces  the  torque  and  speed  requirements  of  the 
motor. 

2.  Relationship  Between  Louver  Angle  and  Open  Area 

The  RWE  open  area  is  determined  by  the  amount  of  louver  area  seen  by  the  expansion 
section.  Imagine  looking  down  the  expansion  section  into  the  RWE.  One  cannot  see  the 
portions  of  the  louver  which  lie  outside  the  circular  area  of  the  expansion  section.  Thus, 
these  portions  of  the  louvers  are  not  considered  when  calculating  the  RWE  open  area.  Only 
that  area  projected  onto  the  louvers  by  the  shock  tube  is  considered. 

The  area  of  the  exit  plane  of  the  RWE  is  always  given  relative  to  the  area  of  the 
expansion  section  of  the  shock  tube.  Therefore,  if  the  exit  plane  area  of  the  RWE  is  equal 
to  the  area  of  the  expansion  section  then  the  RWE  open  area  is  said  to  be  100  percent.  In 
reality  the  open  area  cannot  reach  100  percent  because  the  louvers  have  some  finite  thickness 
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which  is  presented  to  the  expansion  section  even  when  the  louvers  are  horizontal.  This  can 
be  seen  in  Figure  5,  which  is  a  chart  defining  the  RWE  open  area  as  a  function  of  louver 
angle. 

Relationship  Between  Louver  Angle  and  RWE  Open  Area  Ratio 
Method  Used  for  Majority  of  Program 


Louver  Angle  (deg) 


Figure  5.  Relationship  Between  Louver  Angle  and  Open  Area 

The  maximum  open  area  of  the  RWE  is  87  percent.  The  open  area  is  constant  from 
zero  to  about  6  degrees  due  to  the  diamond  shaped  cross  section  of  the  louvers.  A  louver 
must  rotate  approximately  6  degrees  before  its  tips  will  exit  the  area  blocked  by  the  thicker 
mid-section  of  the  louver.  The  flat  portion  of  the  graph  represents  0  percent  open  area,  where 
the  presented  area  of  the  louvers  equals  the  area  of  the  expansion  section  of  the  shock  tube. 
When  the  louvers  begin  to  overlap  (at  about  61  degrees),  the  open  area  goes  to  0  percent 
and  stays  there  until  the  louvers  exit  the  overlap  condition  at  about  124  degrees. 


3.  The  Electric  Servomotor 

The  four  louvers  are  connected  to  a  gear  train  and  driven  by  an  electric  servomotor. 
The  gears  attached  to  the  louvers  are  each  6.35  cm  in  diameter  and  the  gear  attached  to  the 
motor  shaft  is  2.54  cm  in  diameter.  Included  in  the  documentation  that  came  with  the  motor 
was  the  manufacturers  technical  data  about  the  motor  From  this  the  peak  torque,  rated 
speed,  rotor  inertia  and  acceleration  at  peak  torque  were  obtained.  From  the  drawings  of 
the  RWE,  the  inertia  of  the  louvers  and  gears  was  determined.  With  all  of  this  information 
it  was  possible  to  calculate  the  louver’s  maximum  angular  acceleration  and  the  maximum 
angular  velocity.  These  calculations  can  be  found  in  Appendix  A. 
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4.  Triggering  and  Computer  Control 


As  previously  stated,  a  stagnation  probe  mounted  in  the  expansion  section  near  the 
diaphragm  was  used  to  trigger  the  recording  device  and  the  RWE. 

Control  of  the  motor  was  accomplished  through  a  circuit  board  supplied  with  the  motor 
and  placed  into  a  Zenith  248  personal  computer  system.  Communication  with  the  circuit 
board  was  made  possible  via  a  computer  program  (written  in  BASIC)  which  the  end  user 
could  write  to  suit  his  own  needs.  Prior  to  a  shot,  the  louvers  would  be  set  to  a  predetermined 
angle  and  this  program  would  be  run.  The  program  reads  in  a  series  of  execution  commands 
from  an  ASCII  file  and  loads  them  into  the  motor’s  memory.  At  this  point  the  RWE  is 
armed  and  awaiting  the  trigger  signal.  The  driver  is  then  supplied  with  nitrogen  gas  until 
the  diaphragm  is  ruptured.  The  rupturing  of  the  diaphragm  releases  the  driver  gas  and 
initiates  a  shock  in  the  expansion  section.  When  the  shock  impinges  on  the  probe  at  52.1  cm 
it  triggers  the  recording  sequence  and  the  RWE  motor. 


IV.  The  Experimental  Program 


The  purpose  of  an  RWE  is  to  eliminate  the  rarefaction  wave  formed  at  the  exit  plane  of 
a  shock  tube.  If  the  RWE  is  functioning  properly,  the  experimented  records  should  appear 
as  if  the  expansion  section  length  were  infinite.  Therefore  a  series  of  baseline  tests  were 
performed  using  a  very  long  expansion  section  (1719.6  cm).  Using  this  long  expansion 
section,  a  rarefaction  wave  is  still  formed  but  it  does  not  affect  the  recorded  pressure  histories 
until  well  after  the  period  of  interest.  The  data  collected  from  these  long  expansion  section 
tests  is  also  used  to  generate  the  driving  functions  for  the  RWE. 

Another  series  of  tests  was  performed  employing  a  short  expansion  section  (856.6  cm). 
These  tests  were  conducted  to  show  the  effe'  c  of  the  rarefaction  wave  on  the  pressure  histories 
when  no  RWE  is  used. 

For  the  final  test  series  the  RWE  was  attached  to  the  end  of  this  short  expansion  section. 
The  functions  generated  from  the  long  expansion  section  tests  were  fed  into  the  RWE  for 
each  shot.  In  theory,  the  long  expansion  section  test  and  the  RWE  test  should  be  nearly 
identical  for  a  given  driver  condition  while  the  short  expansion  test  with  no  RWE  will  be 
quite  different  due  to  the  effect  of  the  rarefaction  wave. 

In  order  for  the  experimental  program  to  be  successful,  the  set  of  driver  conditions  had 
to  be  nearly  identical  for  each  test  series.  However,  due  to  varying  diaphragm  characteris¬ 
tics  the  driver  conditions  at  diaphragm  rupture  could  not  be  repeated  exactly  from  test  to 
test.  Each  set  of  driver  conditions  was  comprised  of  four  driver  pressure  levels  and  three 
different  driver  volumes,  resulting  in  twelve  tests  per  series.  The  three  driver  volumes  were 
obtained  by  adjusting  the  length  of  the  cylindrical  driver  and  were  2839.9  cm3,  4899.2  cm3 
and  9017.7  cm3.  The  four  driver  pressure  levels  resulted  from  varying  the  thickness  of  the 
diaphragm  which  was  ruptured  naturally.  From  these,  four  different  shock  overpressures 
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were  obtained  and  were  approximately  40  kPa ,  68  kPa,  88  kPa  and  100  kPa. 


V.  Development  of  the  RWE  Driving  Functions 

The  RWE  driving  functions  were  generated  by  a  FORTRAN  computer  program  which 
used  the  data  collected  from  the  long  expansion  section  tests  as  input.  The  objective  of  the 
program  is  to  control  the  motion  of  the  RWE  louvers  in  such  a  way  that  the  blast  wave 
resulting  from  a  test  using  the  RWE  is  nearly  identical  to  the  long  expansion  section  test 
that  generated  that  particular  RWE  function. 

It  should  be  noted  that  this  program  can  also  generate  functions  from  computational 
results.  Currently  the  program  can  employ  results  of  the  BRL-Q1D  code  and  can  be  easily 
modified  to  accept  results  of  other  hydrodynamic  computer  codes  (hydrocodes).  Thus  far, 
the  computational  models  which  have  been  used  have  not  proven  accurate  enough  to  provide 
reliable  input  to  the  RWE  function  generating  program.  For  this  reason,  experimental  data 
was  used  to  generate  the  functions  outlined  in  this  report. 


1.  Theory  of  RWE  Operation 

The  object  of  using  an  RWE  is  the  elimination  or  suppression  of  the  waves  that  are 
generated  when  a  shock  exits  the  downstream  end  of  a  shock  tube.  It  will  not  be  used  to 
deliberately  emit  a  shock  or  rarefaction  wave  in  order  to  tailor  the  flow  in  the  test  section. 
All  tailoring  of  the  flow  in  the  test  section  will  be  accomplished  by  changes  in  the  driver  end 
of  the  shock  tube. 

Using  an  RWE  to  help  modify  a  blast  wave  in  a  shock  tube  involves  some  difficulties. 
First,  the  waves  generated  by  the  RWE  are  not  the  type  of  waves  found  in  free  field  blast. 
The  shocks  or  rarefactions  generated  by  the  RWE  run  upstream  into  the  flow.  In  the  free 
field  case  the  pressure  histories  are  controlled  by  rarefactions  and  secondary  shocks  running 
downstream  and  overtaking  the  primary  shock  which  cause  the  pressure  to  decay  with  time. 
The  influence  of  the  waves  on  the  flow  depends  on  their  direction.  A  rarefaction  overtaking 
the  shock  causes  a  decrease  in  static  pressure,  flow  velocity  and  dynamic  pressure.  An 
upstream  moving  rarefaction  causes  a  decease  in  static  pressure  but  also  causes  an  increase 
in  flow  velocity  and  can  cause  an  increase  in  dynamic  pressure. 

The  second  difficulty  is  that  an  RWE  can  influence  only  part  of  the  flow  history.  Changes 
at  the  driver  end,  on  the  other  hand,  influence  the  entire  history.  The  magnitude  of  the  shock 
or  rarefaction  returned  by  the  RWE  depends  on  the  RWE  open  area  and  the  strength  of  the 
incident  shock.  The  time  at  which  waves  generated  at  the  RWE  arrive  at  the  test  station 
depends  on  the  decay  rate  of  the  shock  wave  as  well  as  the  incident  shock  overpressure  and 
RWE  open  area.  Controlling  the  strength  and  arrival  time  of  the  RWE  generated  wave  to 
produce  a  desired  change  in  the  test  section  flow  is  a  fairly  complicated  task.  Because  of 
these  difficulties  it  was  decided  to  limit  the  goal  of  RWE  operation  to  simply  eliminating 
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waves  generated  at  the  downstream  end  rather  than  using  them  to  control  flow  at  the  test 
station. 

In  order  to  eliminate  the  rarefaction  waves  generated  at  the  downstream  end,  one  of 
two  conditions  must  exist.  Either  the  flow  exiting  the  shock  tube  has  to  have  the  same 
static  pressure  as  the  atmosphere  or  the  flow  exiting  must  be  sonic  with  an  exiting  static 
pressure  equal  to  or  greater  than  the  atmospheric  pressure.  None  of  the  shocks  in  the  design 
envelope  of  the  LB/TS  are  sufficiently  strong  to  produce  sonic  or  supersonic  flow  behind 
them.  Therefore,  the  following  analysis  is  intended  only  for  shock  waves  which  generate 
subsonic  flow. 

The  initial  interaction  of  the  shock  and  RWE  will  not  be  analyzed.  This  initial  inter¬ 
action  produces  a  pressure  spike  that  moves  upstream  and  can  affect  the  flow  in  the  test 
section.  This  effect,  however,  is  small  and  the  addition  it  makes  to  the  positive  phase  impulse 
is  negligible. 

The  following  assumptions  have  been  made  to  simplify  the  analysis: 

1 .  The  RWE  will  be  assumed  to  act  as  a  simple  converging  nozzle. 

2.  The  flow  will  be  considered  one  dimensional. 

3.  The  interaction  between  the  flow  behind  the  shock  and  the  RWE  is  assumed  to  be 

isentropic.  Entropy  is  generated  when  gas  is  processed  by  the  leading  shock  but  the 

flow  behind  the  shock  can  be  reasonably  modeled  as  isentropic. 

4.  The  flow  is  assumed  to  be  inviscid. 

5.  The  gas  involved  is  assumed  to  have  an  ideal  gas,  polytropic,  equation  of  state. 

6.  The  period  of  interest  is  only  that  period  in  which  the  flow  is  exiting  the  shock  tube. 

The  analysis  does  not  consider  flow  from  the  atmosphere  back  into  the  expansion  tunnel. 

For  isentropic  flow  through  a  converging  nozzle  to  reach  sonic  velocity  in  expanding 
from  a  stagnation  pressure  of  Fo«  to  a  static  atmospheric  pressure  of  the  ratio  Poo/Poi 
must  be  larger  than  the  critical  ratio  given  by  Equation  1  5. 


Poo 

Pm 


(1) 


If  the  ratio  of  specific  heats,  7,  is  equal  to  1.400,  this  critical  ratio  is  0.528282.  If  the 
result  of  Equation  1  is  greater  than  0.528282  the  flow  will  remain  subsonic  at  the  exit  of  the 
RWE  converging  nozzle.  If  the  ratio  is  less  than  that,  then  with  the  proper  RWE  setting  the 
flow  will  choke  and  become  sonic  at  the  exit  of  the  RWE  converging  nozzle.  These  are  the 
two  cases  to  be  considered  in  determining  the  open  area  at  the  end  of  the  RWE  converging 
nozzle.  In  one  case  the  flow  remains  subsonic  and  in  the  other  the  flow  becomes  sonic  at  the 
RWE  exit. 
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The  flow  behind  a  72.2  kPa  overpressure  shock  has  a  Poo/Po.  ratio  greater  than  0.528282 
and  will  therefore  have  subsonic  flow  exiting  the  RWE.  For  shock  overpressures  exceeding 
72.2  kPa  the  ratio  is  less  than  0.528282  and  the  flow  exiting  the  RWE  through  the  appropriate 
area  reduction  initially  will  be  sonic. 

The  case  in  which  the  flow  remains  subsonic  will  be  considered  first.  The  Mach  number 
at  the  inlet  to  the  RWE  is  assumed  to  be  the  undisturbed  local  Mach  number  behind  the 
shock.  Using  the  definition  of  stagnation  pressure,  Pq,  and  solving  for  the  flow  Mach  number 
at  the  RWE  inlet,  A/,-  yields  Equation  2  6. 


If  pi  and  Poi  are  the  static  and  stagnation  pressures  respectively  at  the  RWE  inlet  then  the 
Mach  number  at  the  inlet,  Afj,  can  be  calculated  using  Equation  2. 

For  isentropic  flow  the  stagnation  pressure  remains  constant  through  the  converging 
nozzle.  Therefore  the  stagnation  pressure  at  the  RWE  exit,  P0e,  equals  Po,  (the  known 
inlet  stagnation  pressure).  For  the  rarefaction  wave  to  be  eliminated  the  static  pressure  at 
the  RWE  exit  must  match  the  atmospheric  pressure,  p^,  which  is  also  known.  Therefore 
Equation  2  can  be  modified  and  used  to  determine  the  Mach  number  at  the  RWE  exit,  Me 
as  a  function  of  P<*  and  p,*,  as  shown  in  Equation  3. 


The  relationship  between  Mach  number  and  area  ratio  for  an  isentropic  flow  is  given  by 
Equation  4. 

A.  (MT)  [2 +  (7-1  )Ml]  ^  ... 

At  \Mj  [2  +  (7  -  1 ) A/,2 J  [l 

In  Equation  4,  Ai  is  the  known  inlet  area  to  the  RWE.  Ae  is  the  RWE  exit  area  which  will 
produce  a  match  between  the  static  pressure  of  the  exiting  flow  and  the  atmosphere. 

For  a  decaying  blast  wave  simulation,  the  static  and  stagnation  pressures  at  the  inlet 
of  the  RWE  change  with  time.  As  long  as  these  pressure  histories  are  known,  an  RWE  open 
area  history  can  be  found  by  repeated  calculations  using  Equations  2,  3  and  4. 

For  the  second  case,  in  which  the  ratio  Poo/Poi  is  less  than  0.528282,  the  flow  at  the  exit 
plane  of  the  RWE  is  sonic,  thus  Afe  =  1.0.  Equation  2  is  still  used  to  calculate  M,  but  since 
Me  is  known  Equation  3  is  not  used.  In  this  case  Equation  4  simplifies  to  Equation  5. 

4l  =  M  7  +  1  1^ 

Ai  *  [2 +  (7-  1)A/,2 
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(5) 


For  a  decaying  blast  wave  simulation  with  an  initial  Pco/Poi  ratio  less  than  0.528282, 
Equation  5  is  used  from  the  start  of  the  calculation  until  the  value  of  this  ratio  reaches 
0.528282.  When  this  point  is  leached,  Equation  4  is  used  for  all  subsequent  area  ratio 
calculations  in  the  history. 


2.  Adaptation  of  the  RWE  Theory 

The  original  intent  of  the  experimental  program  was  to  use  the  FORTRAN  code  “ELIM”  2 
developed  by  Dr.  J.  Gottlieb  of  the  University  of  Toronto  Space  Institute  to  develop  the  RWE 
functions.  The  ELIM  code  extends  the  theory  given  in  the  preceding  section  to  take  into  ac¬ 
count  edge  effects  and  also  handles  inflow  as  well  as  outflow.  The  BRL-QlD  7  code  was  used 
to  calculate  the  static  pressure,  sound  speed  and  flow  velocity  as  a  function  of  time  for  the 
1:57  scale  blast  simulator  equipped  with  a  very  long  expansion  tunnel  but  no  RWE.  Attempts 
were  then  made  to  use  the  BRL-QlD  code  results  to  produce  RWE  functions  which  would 
control  the  active  RWE  in  experiments  on  the  1:57  scale  blast  simulator.  These  attempts 
proved  unsuccessful.  Comparison  with  actual  experiments  showed  that  while  the  BRL-QlD 
code  provides  a  reasonable  approximation  of  the  flow  conditions,  it  is  not  accurate  enough 
in  its  details  for  proper  RWE  control.  Attempts  using  the  Random  Choice  Method  (RCM) 
code  to  model  the  flow  in  the  1:57  scale  blast  simulator  also  failed  to  produce  acceptable 
functions. 

An  attempt  was  then  made  to  generate  the  RWE  functions  from  experimental  flow  mea¬ 
surements  recorded  during  long  expansion  tunnel  tests.  This  was  recognized  as  an  expedient 
alternative  in  order  to  carry  out  the  small  scale  active  RWE  testing.  In  a  large  scale  blast 
simulator  it  will  be  impossible  to  obtain  long  expansion  tunnel  pressure  histories  experi¬ 
mentally  due  to  the  impracticality  of  changing  tunnel  lengths.  Work  continues  in  finding  a 
hydrocode  with  sufficient  accuracy  to  produce  flow  histories  on  which  RWE  control  functions 
can  be  based. 

The  ELIM  code  was  not  used  to  generate  RWE  functions  from  the  experimental  data. 
The  sound  speed  and  velocity  required  by  ELIM  could  not  be  derived  from  the  experimental 
instrumentation  available  for  this  work.  A  new  FORTRAN  subroutine,  “AREARP”,  was 
developed  based  on  the  theory  outlined  in  the  previous  section  and  requires  only  atmospheric 
conditions  and  the  RWE  inlet  static  and  stagnation  histories  as  input.  The  subroutine 
AREARP  can  use  either  experimental  or  calculated  stagnation  and  static  pressure  histories 
as  input.  The  histories  must  be  for  a  flow  in  which  there  are  no  waves  generated  at  the 
RWE,  (i.e.,  the  histories  that  would  exist  if  the  expansion  tunnel  continued  to  infinity  beyond 
the  RWE  location).  Such  histories  can  be  approximated  by  using  long  expansion  tubes  in 
experiments  or  calculations  which  keep  waves  from  returning  during  the  period  of  interest. 
The  stagnation  and  static  pressure  histories  were  used  to  determine  the  RWE  area,  ratio 
histories  for  subsequent  experiments  which  employed  a  short  expansion  tunnel  and  an  active 
RWE.  Due  to  the  nature  of  the  pressure  histories  and  the  shock  overpressures  obtainable  in 
the  1:57  scale  blast  simulator  these  attempts  also  proved  unsuccessful. 
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In  the  experimental  setup  shock  overpressures  were  limited  to  a  maximum  of  140  kPa. 
Also,  in  order  to  test  the  theory  for  subsonic  flow  through  the  RWE,  part  of  the  testing 
had  to  be  conducted  with  shock  overpressures  less  than  72.2  kPa.  Unfortunately,  at  low 
overpressures  the  ratio  of  p,  to  Pq,  is  very  close  to  unity.  Under  these  conditions  the  random 
noise  which  exists  in  the  experimental  pressure  histories  tends  dominate  the  calculation  of 
inlet  Mach  number,  Mi.  The  resulting  random  variations  in  Mi  in  turn  produce  an  RWE 
area  ratio  history  with  large  random  variation,  which  is  unacceptable. 

R.  Guice  8  of  Applied  Research  Associates,  Inc.,  suggested  that  the  following  equation 
could  be  substituted  in  place  of  Equation  2. 


A/,= 


Pi/Poo  -  1 

iJiii  -  1)/2t)  (p./Poo)  ((7  -  l)/(7+  1)  +  (Pt/Poo)) 


Equation  6  is  normally  used  to  calculate  the  Mach  number  of  the  flow  directly  behind 
a  moving  shock.  For  the  moving  shock  case  p0 0  would  be  the  pressure  ahead  of  the  shock 
and  p,  would  be  the  pressure  behind  the  shock.  Here  we  are  using  the  equation  to  calculate 
Mach  number  not  only  directly  behind  the  shock  but  also  upstream  where  the  pressure  has 
decayed  from  its  initial  value  behind  the  shock  front. 


The  pressure  decay  behind  the  shock  in  the  blast  simulator  is  caused  by  rarefaction 
waves  generated  in  the  driver  which  overtake  the  shock.  The  Mach  number  directly  behind 
the  incident  shock,  M„  is  calculated  using  Equation  7. 


M  = 


PJPoo  ~  1 

iJiil  ~  l)/27)  (pJPoo)  ((7  -  l)/(7  +  1)  +  PjPoo) 


In  Equation  7,  p,  is  the  static  pressure  directly  behind  the  shock.  Mi,  the  local  Mach  number 
of  the  flow  after  its  pressure  has  decayed  from  p,  to  p,,  is  given  by  Equation  8. 


((7  —  l)M,  —  2) 


(8) 


It  has  been  found  that  using  Equation  6  to  calculate  Mi  (after  the  pressure  has  decayed 
behind  the  shock)  provides  a  very  good  estimate.  This  Mach  number  falls  within  1  %  of  the 
value  that  is  calculated  more  formally  by  first  calculating  M ,  using  Equation  7,  then  using 
Equation  8  to  calculate  Mi .  Equation  6  provides  a  good  estimate  only  if  the  flow  behind  the 
shock  is  subsonic  and  the  contact  surface  between  the  driver  gas  and  the  original  expansion 
section  gas  does  not  pass  through  the  measurement  location.  A  relatively  long  distance  was 
used  between  the  diaphragm  and  the  measurement  stations  in  the  test  series  reported  here. 
The  driver  gas  behind  the  contact  surface  therefore  did  not  reach  the  measurement  stations 
during  the  time  of  interest.  A  review  of  the  stagnation  histories  showed  that  the  contact 
surface  passed  though  very  late  in  the  experiments. 

A  new  FORTRAN  subroutine,  “RWE5”,  was  written  based  on  the  use  of  Equation  6  to 
calculate  the  RWE  inlet  Mach  number,  Mt.  Equation  3  uses  the  inlet  stagnation  pressure, 
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Mach  Number  Histories  from  Different  Subroutines 
Static  Overpressure  =  360  kPa 


Figure  6.  Comparison  of  Mach  Number  Histories 


Poi  in  order  to  calculate  RWE  exit  Mach  number,  Me.  Equation  2,  when  solved  for  Poi  in 
terms  of  Mi,  yields  Equation  9. 


Poi  =  Pi 


2  +  (7  —  1)M? 


■1-1 


(9) 


Once  Poi  has  been  calculated  using  Equation  9,  “RWE5”  proceeds  like  “AREARP” .  De¬ 
pending  on  the  ratio  Poo/Poi ,  the  RWE  area  ratio  Ae/Ai  is  calculated  using  either  Equation  4 
or  Equation  5. 

Figure  6  shows  a  comparison  of  the  Mach  numbers  calculated  in  “AREARP”  and 
“RWE5” .  One  can  see  the  results  are  very  close.  Similar  comparisons  for  other  flow  condi¬ 
tions  of  interest  showed  that  the  Mach  numbers  calculated  by  the  two  algorithms  were  also 
very  close. 


3.  The  Overall  RWE  Function 

Several  characteristics  of  the  RWE  motion  profile  were  decided  upon  early  in  the  pro¬ 
gram.  The  most  obvious  of  these  is  the  idea  that  the  final  RWE  function  should  not  attempt 
to  exceed  the  motor’s  torque  and  speed  limitations.  Secondly,  the  direction  of  rotation  of 
the  louvers  should  not  change  during  a  test.  As  a  result  of  these  limitations,  the  final  mo¬ 
tion  profile  fed  to  the  motor  will  not  exactly  match  that  of  the  exact  function  which  was 
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generated  from  the  pressure  data  alone.  The  final  RWE  function  will  adhere  as  closely  as 
possible  to  the  exact  function,  smoothing  through  reversals  in  direction  and  quick  changes 
in  velocity  which  exceed  the  peak  torque  of  the  motor.  The  differences  between  the  exact 
and  final  functions  are  illustrated  in  Figure  7. 

DIFFERENCE  BETWEEN  EXACT  AND  FINAL  RWE  FUNCTIONS 


'  i  i  i - 1 - 1  i 

0.00  0.01  0.02  0.03  0.04  0.05 

Tlm«  Aftar  Shock  Arrival  (s) 


Figure  7.  Exact  and  Final  Open  Area  versus  Time 

It  should  be  noted  that  the  equations  used  to  calculate  the  RWE  open  area  are  only  valid 
during  the  positive  phase  of  the  blast.  Computational  studies  were  performed  to  determine 
the  nature  of  the  open  area  history  after  the  positive  phase.  The  computer  program  ELIM 
developed  by  Dr.  J.  Gottlieb  was  used  in  this  study.  This  study  determined  that  after  the 
positive  phase  the  RWE  open  area  would  oscillate  about  zero  in  a  manner  similar  to  that  of 
the  static  pressure.  A  few  tests  were  performed  using  this  approach.  These  tests  resulted 
in  damaged  louvers.  At  this  point  it  was  decided  that  the  function  should  reopen  the  RWE 
as  quickly  as  possible  following  the  positive  phase.  Therefore,  when  the  louvers  begin  to 
overlap,  the  louvers  are  accelerated  at  the  maximum  rate  until  the  peak  velocity  is  reached. 
This  velocity  is  maintained  until  the  RWE  is  once  again  fully  open.  Then  the  louvers  are 
brought  to  rest  in  the  maximum  open  position. 

Since  the  RWE  is  a  machine  with  torque  limitations,  some  time  interval  is  required 
for  the  motor  to  start  from  rest  and  achieve  the  initial  velocity  defined  in  the  function. 
Therefore,  the  louvers  are  set  in  motion  as  soon  as  the  trigger  signal  senses  the  flow,  and 
the  motor  is  able  to  overcome  the  resistance  to  initial  acceleration  prior  to  the  arrival  of  the 
shock  wave.  The  louvers  are  pre-positioned  in  such  a  way  that  when  the  shock  arrives  at  the 
RWE,  the  louvers  are  at  the  appropriate  angle  to  eliminate  the  rarefaction  wave. 
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4.  Servomotor  Execution  Commands 


The  motion  commands  which  are  sent  to  the  motor  are  simply  a  series  of  angular 
displacement  values  to  be  executed  during  a  fixed  time  interval.  This  angular  displacement 
value  represents  the  angular  distance  to  be  moved  in  that  time  interval,  i.e.,  a  displacement 
relative  to  the  position  at  the  start  of  that  time  interval.  The  commands  are  not  absolute 
commands  to  go  to  some  specific  angular  position  relative  to  a  fixed  origin.  The  fixed  time 
interval  can  be  defined  within  the  RWE  controller  program  discussed  earlier.  For  the  electric 
motor  used,  this  time  interval  must  be  in  the  range  of  2  —  50  ms.  Since  the  blast  waves  in 
question  are  short  duration  events  (less  than  50  ms),  the  minimum  time  interval  of  2  ms 
was  chosen  so  that  the  control  of  the  motor  could  be  as  fine  as  possible. 

The  angular  displacement  value  for  a  particular  time  interval  is  represented  by  the 
amount  of  motor  angular  units  ( MAU )  the  rotor  should  move  in  that  time  period.  It  takes 
25000  MAU  for  one  complete  rotation  (360  degrees)  of  a  louver.  Therefore,  69.4  MAU  are 
required  to  move  the  louvers  1  degree  in  one  time  interval. 

The  series  of  motor  commands  is  generated  by  first  breaking  up  the  ideal  function 
into  even  time  steps  of  2  ms,  interpolating  for  RWE  open  area  where  necessary.  For  each 
area  ratio  in  the  function  there  is  a  correct  louver  angle  which  is  calculated  from  the  RWE 
geometry.  The  change  in  louver  angle  from  one  time  period  to  the  next  is  multiplied  by  69.4 
to  yield  the  number  of  MAU  to  be  executed  for  that  time  interval.  The  calculated  MAU 
value  is  then  rounded  to  the  nearest  integer. 

Calculations  in  Appendix  A  indicate  that  the  motor  should  be  able  to  achieve  an  angular 
velocity  of  46.08  rad/s  and  an  angular  acceleration  of  5687.6  rad/s2.  This  angular  velocity 
limitation  means  that  the  maximum  angular  displacement  that  the  motor  can  achieve  in  a 
two  millisecond  time  interval  is  0.09216  rad  or  366.7  MAU.  To  be  conservative,  the  computer 
code  which  generates  the  motor  commands  limited  the  maximum  angular  displacement  in 
any  2  ms  time  interval  to  342  MAU.  If  the  desired  displacement  exceeded  342  MAU  that 
displacement  value  was  set  to  342  MAU.  The  limit  on  angular  acceleration  means  the 
maximum  change  in  angular  velocity  during  a  2  ms  time  step  is  11.38  rad/s.  The  computer 
code  which  generates  the  commands  was  constructed  in  such  a  way  that  if  the  desired 
acceleration  exceeded  the  maximum,  the  acceleration  was  set  equal  to  the  maximum. 


5.  The  Effect  of  Driver  Volume  on  the  RWE  Function 

The  initial  open  area  of  the  RWE  is  defined  by  the  shock  pressure.  The  positive  phase 
duration  (PPD)  of  the  blast  wave  defines  the  time  required  for  the  RWE  to  change  from 
the  initial  setting  to  the  overlap  (fully  closed)  position.  As  the  driver  volume  is  reduced, 
the  PPD  decreases.  From  this  one  can  see  that  smaller  driver  volumes  require  the  greatest 
angular  velocities  and  accelerations  of  the  RWE. 

Another  problem  posed  by  the  small  driver  volume  is  the  size  of  the  time  increment.  For 
a  blast  wave  whose  PPD  is  50  ms,  25  motor  commands  will  be  executed  during  the  positive 
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phase.  When  the  PPD  is  reduced  to  14  ms,  only  7  commands  will  be  executed  during  the 
positive  phase.  Therefore,  as  the  driver  volume  is  reduced  the  resolution  of  control  becomes 
more  coarse,  possibly  resulting  in  greater  error  between  the  ideal  and  real  functions. 


VI.  Preliminary  Tests 

Before  the  test  series  used  to  validate  the  theory  of  RWE  operation  was  performed,  a 
number  of  preliminary  tests  were  conducted.  Several  tests  were  performed  to  confirm  that 
the  RWE  hardware  and  software  had  been  properly  integrated  with  the  1:57  scale  blast 
simulator.  In  these  early  tests  a  simple  linear  change  in  RWE  area  ratio  was  used.  The 
linear  closing  function  started  with  the  RWE  set  for  the  flow  condition  immediately  behind 
the  shock  and  then  closed  to  0  percent  open  area  linearly  in  a  time  equal  to  the  positive 
phase  duration  for  the  blast  wave  in  the  simulator.  The  linear  change  in  RWE  open  area 
produced  acceptable  results  in  the  few  preliminary  tests  in  which  it  was  used.  The  object 
of  this  study,  however,  was  to  validate  a  general  theory  of  RWE  operation,  not  to  obtain 
empirically  derived  RWE  closing  functions.  For  this  reason  linear  closing  functions  were  not 
pursued  further. 

The  next  set  of  tests  attempted  to  use  flow  field  data  generated  by  the  BRL-QlD  code 
as  input  to  the  program  which  generates  the  RWE  closing  functions.  The  BRL-QlD  code 
was  used  to  simulate  the  flow  in  the  1:57  scale  blast  simulator  equipped  with  a  very  long 
expansion  section.  The  purpose  of  this  exercise  was  to  eliminate  the  necessity  of  actually 
performing  tests  in  the  1:57  scale  blast  simulator  with  an  additional  expansion  tunnel  length 
downstream  of  the  normal  RWE  location.  The  flow  field  data  generated  by  the  BRL-QlD 
code,  however,  was  not  accurate  enough  to  produce  acceptable  functions.  After  several  un¬ 
successful  attempts,  this  method  was  abandoned.  In  all  subsequent  tests,  the  RWE  functions 
were  generated  from  data  gathered  in  experiments  performed  with  an  additional  expansion 
tunnel  length  added  downstream  of  the  normal  RWE  location. 

Many  tests  were  performed  to  determine  the  best  way  of  modeling  the  geometry  of  the 
flow  passing  through  the  RWE.  At  one  extreme  the  flow  can  be  thought  of  as  filling  the  entire 
cross-section  of  the  RWE  which  has  a  rectangular  cross-sectional  area  larger  than  that  of 
the  circular  expansion  section.  In  the  other  case,  the  flow  is  considered  to  form  a  jet  which 
retains  the  shape  of  the  original  expansion  tunnel.  The  real  case  actually  lies  somewhere 
between  the  two  extremes  and  changes  with  time  during  a  test.  In  the  preliminary  tests,  it 
was  assumed  that  the  flow  filled  the  entire  cross-section  of  the  RWE.  The  area  of  the  frame 
and  the  entire  lengths  of  the  louvers  were  considered  when  calculating  the  open  area  of  the 
RWE.  The  relationship  between  louver  angle  and  open  area  for  this  method  is  shown  in 
Figure  8.  Since  the  area  of  the  RWE  frame  is  greater  than  that  of  the  expansion  section ,  the 
maximum  open  area  for  this  method  is  greater  than  100  percent.  Also,  the  tolerance  between 
the  louvers  and  the  frame  causes  the  minimum  open  area  to  be  greater  than  0  percent. 

Several  tests  were  performed  using  RWE  closing  functions  derived  on  the  assumption 
that  the  flow  filled  the  entire  cross-section  of  the  RWE.  A  pressure  history  from  one  such 
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Relationship  Between  Louver  Angle  and  RWE  Open  Area  Ratio 
Method  Used  Early  in  Program 


Figure  8.  The  Relationship  between  Louver  Angle  and  Open  Area  (Old) 

test  is  compared  to  a  long  expansion  section  test  in  Figure  9.  The  pressure  histories  in  this 
figure  were  recorded  at  the  station  located  7.519  m  from  the  diaphragm.  In  this  figure,  one 
can  see  that  the  greater  louver  angle  associated  with  this  method  generated  an  upstream 
travelling  compression  wave  which  caused  the  flow  to  back  up  in  the  shock  tube.  Also,  the 
louvers  in  the  RWE  were  damaged  by  the  flow  during  some  of  these  tests.  From  these  results 
it  is  obvious  that  this  was  not  the  correct  method  of  calculating  the  open  area  of  the  RWE. 
The  better  method  is  that  which  assumes  the  flow  retains  the  shape  of  the  expansion  section. 
This  alternative  method,  discussed  earlier,  was  used  in  all  subsequent  tests. 

In  the  preliminary  testing  it  was  found  that  the  RWE  was  incapable  of  following  the  ideal 
function  through  the  negative  pressure  phase  of  the  flow.  In  light  of  this  it  was  decided  to 
decelerate  the  louvers  after  they  had  overlapped  and  bring  them  to  rest  to  the  closed  position. 
This  procedure  was  intended  to  reduce  “organ  piping”  (repeated  passage  of  secondary  waves 
up  and  down  the  simulator).  Holding  the  RWE  closed,  however,  resulted  in  damaged  louvers. 
In  all  tests  to  follow  the  louvers  were  rotated  rapidly  through  the  overlap  condition  and 
brought  to  rest  in  the  maximum  open  position. 

Finally,  a  series  of  tests  was  conducted  in  which  the  RWE  functions  used  were  derived 
by  applying  a  curve  fit  to  the  data  gathered  from  long  expansion  section  tests.  This  was 
an  attempt  to  reduce  the  torque  requirements  placed  on  the  motor  by  candidate  functions. 
Unfortunately,  the  RWE  closing  functions  based  on  this  method  were  not  satisfactory. 

In  all,  the  preliminary  tests  far  outnumbered  the  tests  performed  in  the  main  test 
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Test  Result  Using  Old  Method  of 
Computing  RWE  Open  Areo  Ratio 


Figure  9.  Pressure  History  from  Old  Louver  Angle/Open  Area  Relation 

series.  The  preliminary  tests  played  a  valuable  role  in  developing  a  workable,  if  not  optimal, 
methodology  for  operating  the  1:57  scale  RWE.  The  objective  of  this  study  was  to  validate 
a  general  theory  of  RWE  operation  for  use  in  designing  other  blast  simulators,  particularly 
the  proposed  LB/TS  facility.  The  following  results  section  will  concentrate  on  the  main  test 
series  which  grew  out  of  the  knowledge  gained  in  the  preliminary  tests  and  is  the  vehicle  by 
which  the  general  theory  was  validated. 


VII.  Discussion  of  the  Results 


Pressure-time  histories  are  the  primary  output  collected  from  the  tests.  The  effective¬ 
ness  of  the  RWE  can  be  determined  by  overlaying  the  pressure-time  histories  from  a  long 
expansion  section  test,  a  short  expansion  section  test  and  an  RWE  test  for  a  given  set  of 
driver  conditions.  Positive  phase  duration  and  impulse  can  also  be  determined  from  these 
pressure-time  histories.  Comparison  of  these  values  is  another  way  of  determining  the  effec¬ 
tiveness  of  the  RWE. 

The  main  test  series  consisted  of  12  different  driver  conditions  and  3  expansion  section 
configurations  resulting  in  a  total  of  36  shots.  The  shock  tube  configuration  and  initial 
conditions  for  these  shots  are  outlined  in  Table  1. 

As  stated  earlier,  the  measurement  station  located  7.519  m  from  the  diaphragm  is 
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important  because  the  scaled  distance  from  this  station  to  the  RWE  is  equivalent  to  that 
in  the  full  scale  LB/TS.  Therefore,  the  data  collected  at  this  station  ultimately  determines 
the  success  or  failure  of  the  RWE  in  eliminating  the  rarefaction  waves.  For  this  reason,  the 
emphasis  of  this  discussion  will  be  placed  on  data  collected  at  this  station.  A  complete  set 
of  pressure-time  histories  from  all  of  the  stations  is  presented  in  Appendix  B. 

The  ability  of  the  RWE  to  eliminate  the  rarefaction  wave  is  diminished  as  driver  volume 
is  reduced.  For  this  reason  the  experimental  results  are  discussed  in  order  of  decreasing  driver 
volume. 


1.  Results  for  the  9018  cm3  Driver 

The  results  of  tests  using  driver  a  volume  of  9018  cm3  are  shown  in  Tables  2,  3  and 
4.  The  positive  phase  duration  for  long  expansion  section  shots  using  this  driver  averaged 
about  50  ms.  As  a  result,  there  are  about  25  data  points  in  the  RWE  function  during  the 
positive  phase.  Also,  the  long  positive  phase  duration  means  the  average  angular  velocity 
will  not  be  as  great  as  that  of  a  smaller  driver.  As  a  result,  the  RWE  should  be  able  to 
follow  the  input  function  more  closely. 

Figures  10  through  17  are  graphs  showing  the  results  of  shots  using  the  9018  cm3 
driver.  Each  pair  of  graphs  corresponds  to  one  of  the  four  pressure  levels.  Each  pressure 
level  corresponds  to  a  set  of  tests  performed  using  a  particular  thickness  of  diaphragm. 
Pressure  level  1  corresponds  to  the  thinnest  diaphragm  while  pressure  level  4  represents  the 
thickest.  The  first  graph  in  each  set  compares  the  input  RWE  function  to  the  measured 
RWE  motion  for  that  pressure  level.  The  second  graph  in  a  set  overlays  pressure  histories 
from  long  expansion,  short  expansion  and  active  RWE  tests.  In  all  the  graphs,  the  time 
scale  is  given  relative  to  the  time  of  shock  arrival  at  the  station  located  8.281  m  from  the 
diaphragm. 

Overall,  the  RWE  was  able  to  follow  the  input  function  for  the  four  cases.  However,  for 
the  tests  at  pressure  levels  2,  3  and  4  the  RWE  seems  to  overshoot  the  input  function  late 
in  the  positive  phase.  For  each  of  these  cases  the  area  ratio  history  appears  to  be  correct 
during  the  early  part  of  the  positive  phase.  The  greatest  deviation  in  the  area  ratio  occurs 
late  in  the  positive  phase.  In  spite  of  this  difference,  the  resulting  pressure-time  histories 
show  excellent  agreement  between  the  long  expansion  section  shots  and  the  active  RWE 
shots. 

The  effect  of  the  rarefaction  wave  can  be  seen  in  the  short  expansion  section  traces. 
While  this  trace  initially  follows  that  of  the  long  expansion  section,  the  rarefaction  wave 
eventually  passes  through  the  test  station  causing  the  two  traces  to  deviate.  When  this 
happens  the  overpressure  quickly  goes  negative,  then  slowly  returns  to  zero.  As  a  result, 
the  PPD  and  the  impulse  for  a  short  expansion  section  test  are  both  much  less  than  those 
associated  with  a  long  expansion  section  test. 

Residing  in  the  active  RWE  pressure-time  histories  are  several  short  duration,  positive 
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pressure  spikes.  These  transient  spikes  are  caused  by  reflections  from  the  louvers  of  the 
RWE.  This  phenomenon  is  inherent  to  any  rarefaction  wave  eliminator  and  can  never  be 
completely  eliminated.  Thus,  the  resulting  impulse  of  the  RWE  test  is  slightly  greater  than 
the  impulse  of  the  baseline  test. 

Figures  18  and  19  compare  the  impulse  and  PPD  respectively  for  all  the  tests  using 
the  9018  cm3  driver.  In  the  first  of  these  two  charts,  the  long  expansion  section  impulse 
and  the  active  RWE  impulse  are  similar  with  the  latter  slightly  higher  than  the  former,  as 
one  would  expect.  The  positive  phase  durations  for  these  cases  compare  just  as  well  in  the 
bottom  chart.  In  fact  they  are  exactly  the  same  for  the  highest  pressure.  As  expected,  the 
impulse  and  the  PPD  data  for  the  short  expansion  section  tests  fall  well  below  the  results  of 
the  other  cases. 
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Table  1.  Test  Matrix 


Test 

Driver 

Expansion 

RWE 

Ambient 

Ambient 

Driver 

Name 

Volume 

Length 

Used 

Pressure 

Temperature 

Overpressure 

(cm3) 

(cm) 

(-) 

( kPa ) 

(C) 

(. kPa ) 

RWE92 

2839.9 

1719.6 

mmm 

13.9 

1276 

RWE93 

2839.9 

1719.6 

MM 

13.9 

2344 

RWE94 

2839.9 

1719.6 

m 

13.9 

3241 

RWE95 

2839.9 

1719.6 

IBfl 

13.9 

3930 

RWE100 

2839.9 

856.6 

No 

102.9 

17.8 

1207 

RWE101 

2839.9 

856.6 

No 

102.9 

17.8 

2241 

RWE102 

2839.9 

856.6 

No 

102.9 

17.8 

3172 

RWE103 

2839.9 

856.6 

No 

102.9 

17.8 

3792 

RWE96 

2839.9 

856.6 

Yes 

102.9 

17.8 

1276 

RWE97 

2839.9 

856.6 

Yes 

102.9 

17.8 

2517 

RWE98 

2839.9 

856.6 

Yes 

102.9 

17.8 

3172 

RWE99 

2839.9 

856.6 

Yes 

102.9 

17.8 

4068 

RWE61 

4988.2 

1719.6 

No 

mmwm 

1069 

RWE63 

4988.2 

1719.6 

No 

HI 

2310 

RWE65 

4988.2 

1719.6 

No 

MM 

18.3 

3034 

RWE67 

4988.2 

1719.6 

No 

HI 

18.3 

4137 

RWE68 

4988.2 

856.6 

No 

HI 

18.3 

1138 

RWE71 

4988.2 

856.6 

No 

104.1 

18.3 

2344 

RWE72 

4988.2 

856.6 

No 

104.1 

18.3 

3068 

RWE74 

4988.2 

856.6 

No 

104.1 

18.3 

3965 

RWE77 

4988.2 

856.6 

Yes 

104.7 

15.0 

1172 

RWE79 

4988.2 

856.6 

Yes 

104.7 

15.0 

2310 

RWE80 

4988.2 

856.6 

Yes 

104.7 

15.0 

3137 

RWE82 

4988.2 

856.6 

Yes 

104.7 

15.0 

3758 

RWEBL2 

1719.6 

No 

103.5 

22.2 

1482 

RWEBL3 

1719.6 

No 

103.5 

22.2 

2310 

RWEBL5 

9017.7 

1719.6 

No 

103.5 

22.2 

3172 

RWEBL7 

9017.7 

1719.6 

No 

103.5 

22.2 

4275 

RWES2 

9017.7 

856.6 

No 

103.5 

23.9 

1586 

RWES3 

9017.7 

856.6 

No 

103.5 

23.9 

2413 

RWES5 

9017.7 

856.6 

No 

103.5 

23.9 

3275 

RWES7 

9017.7 

856.6 

No 

103.5 

23.9 

4206 

RWE85 

9017.7 

856.6 

Yes 

102.5 

21.7 

1379 

RWE87 

9017.7 

856.6 

Yes 

102.5 

21.7 

2551 

RWE88 

9017.7 

856.6 

Yes 

102.5 

21.7 

3241 

RWE90 

9017.7 

856.6 

Yes 

102.5 

21.7 

3965 
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Table  2.  Static  Overpressure  for  Driver  Volume  =  9018  cm3 


Diaphragm 

Thickness 

M 

Static  Overpressure 
( kPa ) 

Station  Location  =  7.519  m  from  Diaphragm 

Long 

Expansion 

Short 

Expansion 

RWE 

Average 

0.010 

44.349 

47.751 

45.278 

45.793 

0.020 

71.431 

70.609 

69.380 

70.473 

0.030 

89.676 

87.607 

88.984 

88.756 

0.040 

107.724 

107.157 

105.895 

106.925 

Table  3.  Static  Overpressure  Impulse  for  Driver  Volume  =  9018  cm3 


Average 

Static 

Overpressure 

(kPa) 

Static  Overpressure  Impulse 
( kPa  —  s) 

Station  Location  =  7.519  m  from  Diaphragm 

Long 

Expansion 

Short 

Expansion 

RWE 

45.793 

0.515 

0.262 

0.585 

70.473 

0.778 

0.456 

0.894 

88.756 

1.099 

0.628 

1.255 

106.925 

1.348 

0.810 

1.422 

Table  4.  Positive  Phase  Duration  for  Driver  Volume  =  9018  cm3 


Average 

Static 

Overpressure 
(i kPa ) 

Positive  Phase  Duration 
(ms) 

Station  Location  =  7.519  m  from  Diaphragm 

Long 

Expansion 

Short 

Expansion 

RWE 

45.793 

MSM 

9.0 

48.0 

70.473 

9.0 

41.0 

88.756 

54.0 

12.0 

47.0 

106.925 

46.0 

12.0 

46.0  | 
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RWE  Open  Area  Rotio  History 
Driver  Volume  =  9018  cm*  :  Pressure  Level  1 


Driver  Volume  =  9018  cm*  :  Pressure  Level  1 
Test  Station  Location  =  7.519  meters  from  Diaphragm 
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RWE  Open  Area  Ratio  History 
Driver  Volume  =  9018  cm1  :  Pressure  Level  2 


Figure  12.  RWE  Area  History  for  9018  cm3  Driver;  Pressure  Level  2 


Driver  Volume  =  9018  cm1  :  Pressure  Level  2 
Test  Station  Location  =  7.519  meters  from  Diaphragm 
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Overpressure  (kPa) 


RWE  Open  Area  Ratio  History 
Driver  Volume  =  9018  cm*  :  Pressure  Level  3 


Figure  14.  RWE  Area  History  for  9018  cm3  Driver;  Pressure  Level  3 


Driver  Volume  *  9018  cm*  :  Pressure  Level  3 
Test  Station  Location  =  7.519  meters  from  Diaphragm 


25 


Overpressure  (kPa)  QTQ  RWE  Open  Area  RaHo 


RWE  Open  Area  Ratio  History 
Driver  Volume  =  9018  cm’  :  Pressure  Level  4 


-.02  0.00  0.02  0.04  0.06  0.08  0.10  0.12 

Time  After  Shock  Arrival  (s) 


are  16.  RWE  Area 


History  for  9018  cm3 


Driver;  Pressure  Level  4 


Driver  Volume  =  9018  cm*  :  Pressure  Level  4 


Time  After  Shock  Arrival  (s) 

Figure  17.  Pressure  Histories  for  9018  cm3  Driver;  Pressure  Level  4 
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Static  Overpressure  Impulse  (kPa-s) 


Results  of  Experiments  Using  Driver  Volume  =  9018  cm* 
Test  Station  Location  =  7.519  meters  from  Diaphragm 


Figure  18.  Comparison  of  Impulse  for  9018  cm3  Driver 


Results  of  Experiments  Using  Driver  Volume  =  9018  cm* 
Test  Station  Location  *  7.519  meters  from  Diaphragm 


Driver  Pressure  Level 


Figure  19.  Comparison  of  Positive  Phase  Duration  for  9018  cm3  Driver 
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2.  Results  for  the  4988  cm3  Driver 


The  results  of  experiments  using  the  4988  cm3  driver  are  presented  in  the  same  format 
as  those  of  the  previous  driver.  The  experimental  results  for  this  driver  are  shown  in  Tables  5. 
6  and  7.  Figures  20  through  27  show  the  time  histories  obtained  from  some  of  the  tests. 
The  average  long  expansion  section  positive  phase  duration  was  about  30  ms.  Therefore, 
the  demand  placed  on  the  RWE  by  the  functions  should  be  greater  for  this  driver  than  the 
larger  one. 

For  this  driver  configuration  the  RWE  followed  the  function  well  with  the  greatest  error 
occurring  late  in  the  positive  phase.  The  pressure-time  histories  also  compare  well  for  the 
long  expansion  and  the  RWE  tests.  Figures  28  and  29  illustrate  the  comparisons  of  impulse 
and  positive  phase  duration.  As  expected,  the  impulse  values  for  the  long  expansion  and 
RWE  tests  were  similar  but  not  as  close  as  those  of  the  9018  cm3  driver,  and  the  PPD  for 
the  short  expansion  tests  were  well  below  the  others. 

It  is  interesting  to  note  that  as  the  driver  volume  decreases,  the  difference  between  the 
short  expansion  and  long  expansion  impulse  values  also  decreases.  This  occurs  because  the 
flow  in  the  test  section  has  a  shorter  duration  due  to  the  decreased  driver  volume,  leaving 
less  time  for  the  wave  returning  from  the  RWE  to  influence  it. 
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Table  5.  Static  Overpressure  for  Driver  Volume  =  4988  cm3 


Diaphragm 

Thickness 

M 

Static  Overpressure 
(kPa) 

Station  Location  =  7.519  in  from  Diaphragm 

Long 

Expansion 

Short 

Expansion 

RWE 

Average 

■ 

31.905 

33.785 

31.821 

32.504 

68.293 

68.157 

68.180 

68.210 

83.917 

88.269 

85.159 

85.782 

1  0.040 

92.712 

95.898 

96.091 

94.900 

Table  6.  Static  Overpressure  Impulse  for  Driver  Volume  =  4988  cm3 


Average 

Static 

Overpressure 

(kPa) 

Static  Overpressure  Impulse 
(kPa  —  s) 

Station  Location  =  7.519  m  from  Diaphragm 

Long 

Expansion 

Short 

Expansion 

RWE 

32.504 

■rare— 

MMm 

0.273 

68.210 

m 

mm 

0.571 

|  85.782 

■IB 

in 

0.763  | 

94.900 

0.731 

■BB 

0.811 

Table  7.  Positive  Phase  Duration  for  Driver  Volume  =  4988  cm3 


Average 

Static 

Overpressure 

(kPa) 

Positive  Phase  Duration 
(ms) 

Station  Location  =  7.519  m  from  Diaphragm 

Long 

Expansion 

Short 

Expansion 

RWE 

32.504 

■KSB 

30.0 

68.210 

| 

H$B 

30.0 

85.782 

33.0 

36.0 

94.900 

29.0 

9.0 

37.0 
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Static  Overpressure  (kPa) 


RWE  Open  Area  Ratio  History 
Driver  Volume  =  4988  cm*  :  Pressure  Level  1 


Figure  20.  RWE  Area  History  for  4988  cm3  Driver;  Pressure  Level  1 


Driver  Volume  =  4988  cm*  :  Pressure  Level  1 
Test  Station  Location  =  7.519  meters  from  Diaphragm 


Figure  21.  Pressure  Histories  for  4988  cm3  Driver;  Pressure  Level  1 
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RWE  Open  Area  Ratio  History 
Driver  Volume  =  4988  cm1  :  Pressure  Level  2 


Figure  22.  RWE  Area  History  for  4988  cm3  Driver;  Pressure  Level  2 


Driver  Volume  ss  4988  cm*  :  Pressure  Level  2 
Test  Station  Location  =  7.519  meters  from  Diaphragm 


Figure  23.  Pressure  Histories  for  4988  cm3  Driver;  Pressure  Level  2 
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Overpressure  (kPa) 


RWE  Open  Area  Ratio  History 
Driver  Volume  =  4988  cm’  :  Pressure  Level  3 


Driver  Volume  =  4988  cm*  :  Pressure  Level  3 
Test  Station  Location  =  7.519  meters  from  Diaphragm 
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Overpressure  (kPa) 


RWE  Open  Area  Ratio  History 
Driver  Volume  =  4988  cm*  :  Pressure  Level  4 


Figure  26.  RWE  Area  History  for  4988  cm3  Driver;  Pressure  Level  4 


Driver  Volume  a  4988  cm*  :  Pressure  Level  4 
Test  Station  Location  *  7.519  meters  from  Diaphragm 


Figure  27.  Pressure  Histories  for  4988  cm3  Driver;  Pressure  Level  4 
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Static  Overpressure  Impulse  (IcPa-s) 


Results  of  Experiments  Using  Driver  Volume  =  4988  cm1 
Test  Station  Location  =  7.519  meters  from  Diaphragm 


Driver  Pressure  Level 


Figure  28.  Comparison  of  Impulse  for  4988  cm3  Driver 


Results  of  Experiments  Using  Driver  Volume  s  4988  cm* 
Test  Station  Location  =  7.519  meters  from  Diaphragm 


Driver  Pressure  Level 

Figure  29.  Comparison  of  Positive  Phase  Duration  for  4988  cm3  Driver 
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3.  Results  for  the  2840  cm3  Driver 


Finally,  the  results  of  the  tests  using  the  2840  cm3  driver  are  presented  in  Tables  8,  9 
and  10  and  in  Figures  30  through  37.  Again,  the  format  is  the  same  as  the  previous  sets. 
The  PPD  for  the  long  expansion  section  tests  using  this  driver  averaged  15  ms.  This  short 
PPD  indicated  that  the  RWE  would  function  only  marginally  as  intended.  This  supposition 
was  proved  correct  in  the  experiments  using  the  active  RWE. 

The  charts  in  Figures  38  and  39  show  the  marginal  success  of  the  RWE  in  these  tests. 
The  impulse  values  of  the  RWE  tests  compare  well  with  those  of  the  long  expansion  section 
tests.  However,  the  positive  phase  duration  results  are  quite  different.  This  means  the  RWE 
could  not  close  rapidly  enough  to  eliminate  the  rarefaction  wave  for  two  of  the  four  pressure 
levels. 

The  pressure-time  histories  in  Figures  33  and  37  show  a  slight  rarefaction  formed  in  the 
RWE  trace  of  each  figure.  It  is  also  interesting  to  note  that  for  these  shots,  the  greatest  error 
in  the  RWE  history  occurs  in  the  early  part  of  the  positive  phase  as  illustrated  in  Figure  32. 
In  most  of  the  other  active  RWE  tests,  the  RWE  history  was  accurate  early  in  the  positive 
phase  but  deviated  from  the  input  function  late  in  the  positive  phase.  Comparison  of  these 
figures  indicates  that  in  order  to  eliminate  a  rarefaction  wave  the  accuracy  of  the  RWE  area 
is  most  critical  early  in  the  positive  phase.  Louver  position  errors  which  occur  late  in  the 
positive  phase  have  little  effect  on  the  resulting  pressure  history. 
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Table  8.  Static  Overpressure  for  Driver  Volume  =  2840  cm3 


Diaphragm 

Thickness 

(*'») 

Static  Overpressure 
( kPa ) 

Station  Location  =  7.519  m  from  Diaphragm 

Long 

Expansion 

Short 

Expansion 

RWE 

Average 

■SSI 

34.865 

42.074 

42.651 

El 

67.271 

70.515 

66.277 

RUB 

74.481 

74.264 

76.871 

1  0.040 

101.025 

86.395 

89.187 

92.202 

Table  9.  Static  Overpressure  Impulse  for  Driver  Volume  =  2840  cm3 


Average 

Static 

Overpressure 

(kPa) 

Static  Overpressure  Impulse 
( kPa  —  s) 

Station  Location  =  7.519  m  from  Diaphragm 

Long 

Expansion 

Short 

Expansion 

RWE 

42.651 

msSM 

mssm 

0.179 

66.277 

B9 

WmM 

0.257 

76.871 

Bff 

■9 

0.441 

|  92.202 

0.457 

— 

0.444 

Table  10.  Positive  Phase  Duration  for  Driver  Volume  =  2840  cm3 


Average 

Static 

Overpressure 
(i kPa ) 

Positive  Phase  Duration 
(ms) 

Station  Location  =  7.519  m  from  Diaphragm 

Long 

Expansion 

Short 

Expansion 

RWE 

42.651 

15.0 

7.5 

18.0 

66.277 

16.5 

7.5 

9.0 

76.871 

15.5 

8.5 

22.5 

92.202 

13.5 

9.0 

10.0 
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RWE  Open  Areo  Ratio  History 
Driver  Volume  =  2840  cm1  :  Pressure  Level  1 


Driver  Volume  =  2840  cm*  :  Pressure  Level  1 
Test  Station  Location  *  7.519  meters  from  Diaphragm 


Figure  31.  Pressure  Histories  for  2840  cm3  Driver;  Pressure  Level  1 
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Overpressure  (kPa) 


RWE  Open  Area  Ratio  History 
Driver  Volume  =  2840  cm1  :  Pressure  Level  2 


Figure  32.  RWE  Area  History  for  2840  cm3  Driver;  Pressure  Level  2 


Driver  Volume  =  2840  cm*  :  Pressure  Level  2 
Test  Station  Location  a  7.519  meters  from  Diaphragm 


Time  After  Shock  Arrival  (s) 


Figure  33.  Pressure  Histories  for  2840  cm3  Driver;  Pressure  Level  2 


38 


Overpressure  (kPa) 


RWE  Open  Area  Ratio  History 
Driver  Volume  =  2840  cm’  :  Pressure  Level  3 


Figure  34.  RWE  Area  History  for  2840  cm3  Driver;  Pressure  Level  3 


Driver  Volume  =  2840  cm*  :  Pressure  Level  3 
Test  Station  Location  -  7.519  meters  from  Diaphragm 
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RWE  Open  Area  Ratio  History 
Driver  Volume  =  2840  cm*  :  Pressure  Level  4 


Figure  36.  RWE  Area  History  for  2840  cm3  Driver;  Pressure  Level  4 


0 river  Volume  =s  2840  cm*  :  Pressure  Level  4 
Test  Station  Location  =  7.519  meters  from  Diaphragm 


Figure  37.  Pressure  Histories  for  2840  cm3  Driver;  Pressure  Level  4 
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Static  Overpressure  Impulse  (kPa-s) 


Results  of  Experiments  Using  Driver  Volume  =  2840  cm* 
Test  Station  Location  =  7.519  meters  from  Diaphragm 


Driver  Pressure  Level 


Figure  38.  Comparison  of  Impulse  for  2840  cm3  Driver 

Results  of  Experiments  Using  Driver  Volume  =  2840  cm* 

Test  Station  Location  =  7.519  meters  from  Diaphragm 


£  20 


ffn  SHORT 

rrm  long 

res?!  RWE 
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a  5- 
</> 


Driver  Pressure  Level 


Figure  39.  Comparison  of  Positive  Phase  Duration  for  2840  cm3  Dri/er 
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4.  Improving  on  the  Results 


The  experimental  results  showed  that  in  many  tests  the  RWE  failed  to  follow  the  input 
function  at  certain  times  during  the  test.  In  many  cases  this  tracking  error  had  little  effect 
on  the  resulting  pressure-time  histories,  while  in  some  cases  it  had  a  significant  effect.  These 
tracking  errors  can  be  attributed  to  three  causes: 

1.  The  motor  was  not  able  to  accelerate  the  louvers  quickly  enough  to  follow  the  input 
function. 

2.  The  motor  accelerated  the  louvers  properly  but  was  unable  to  decelerate  enough  to 
follow  the  input  function,  resulting  in  the  overshoot  discussed  with  the  results  of  the 
tests  with  the  9018  cm3  driver. 

3.  The  2  ms  sampling  rate  of  the  motor  was  too  coarse  to  provide  accurate  following  of 
the  input  function. 

The  first  two  of  these  three  causes  could  be  overcome  by  using  a  similar  motor  but  one 
with  greater  torque  and  speed  capabilities.  The  solution  to  the  third  cause  is  to  find  a  motor 
with  a  better  sampling  rate  (0.25  —  0.50  ms).  The  important  point  here  is  that  all  of  the 
problems  experienced  with  the  RWE  were  problems  that  could  have  been  solved  by  changing 
or  improving  the  equipment.  No  physically  insurmountable  problems  were  discovered  during 
the  course  of  the  testing.  That  was  the  overall  goal  of  this  project;  not  to  build  a  perfect 
small  scale  RWE,  but  to  prove  that  the  flow  theory  behind  the  development  of  this  RWE 
was  sound,  and  also  to  provide  a  knowledge  base  for  larger  scale  RWE  work.  In  this  respect 
the  program  is  a  success. 


VIII.  Conclusions 

Experiments  have  been  performed  to  determine  the  effectiveness  of  an  active  rarefaction 
wave  eliminator  (RWE)  in  controlling  the  gas  flow  from  the  exit  of  a  shock  tube.  The  purpose 
of  the  tests  was  to  validate  previously  developed  RWE  flow  theory.  The  shock  tube  used  was 
the  1:57  scale  blast  simulator  intended  to  study  the  flow  phenomena  that  will  occur  in  the 
Large  Blast/Thermal  Simulator  being  developed  by  the  Defense  Nuclear  Agency.  Use  of  a 
rarefaction  wave  eliminator  is  necessary  in  order  to  simulate  high  yield  nuclear  events  with 
a  relatively  short  shock  tube. 

The  tests  outlined  in  this  report  employed  three  different  driver  volumes,  four  driver 
pressure  levels  and  three  expansion  section  configurations,  resulting  in  a  total  of  36  shots. 
Success  or  failure  of  the  RWE  was  determined  by  comparing  results  of  active  RWE  tests 
with  similar  long  expansion  section  (no  RWE)  tests. 

The  test  results  show  that  the  RWE  did  indeed  eliminate  the  rarefaction  wave  formed 
at  the  exit  plane  of  the  shock  tube.  The  results  obtained  using  the  two  largest  driver  volumes 
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showed  good  agreement  between  long  expansion  section  section  shots  and  shots  using  the 
RWE.  The  results  gathered  from  the  smallest  driver  volume  tests  yielded  marginal  results. 
The  short  positive  phase  duration  associated  with  this  volume  had  the  RWE  operating  at 
or  beyond  its  limits  of  torque  and  velocity. 

The  final  RWE  functions  outlined  in  this  report  were  the  result  of  many  iterations  of 
test,  analyze  and  re-test.  Numerous  preliminary  tests  were  conducted  before  the  main  test 
series  results  were  obtained.  In  particular,  a  great  deal  of  experimentation  was  performed 
before  the  correct  relationship  between  louver  angle  and  RWE  open  area  ratio  could  be 
realized. 

The  scale  on  which  this  program  was  performed  allowed  for  relatively  simple  changes 
in  shock  tube  configuration  and  rapid  turnaround  between  tests.  This  allowed  us  to  learn  a 
great  deal  in  a  very  short  time.  This  new  knowledge  will  be  a  valuable  asset  when  working 
on  the  1:6  scale  and  full  scale  rarefaction  wave  eliminators. 

The  theory  which  was  validated  by  these  experiments  is  not  influenced  by  the  size  of  the 
system  being  tested.  Therefore,  this  same  theory  should  work  for  the  1:6  scale  and  full  scale 
systems.  However,  parameters  such  as  bearing  friction,  louver  inertia  and  material  stresses 
will  not  scale.  It  is  parameters  such  as  these  that  drive  the  design  of  the  control  system  for 
these  devices.  Simply  scaling  up  the  size  of  the  electric  servomotor  used  in  these  experiments 
by  a  factor  of  57  will  not  yield  an  appropriate  control  system  for  the  full  scale  RWE. 
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APPENDIX  A:  Velocity  and  Acceleration  Limit  Calculations 
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Before  the  closing/opening  functions  for  the  RWE  could  be  derived,  the  capabilities 
of  the  machine  had  to  be  determined.  Once  the  functional  limits  of  the  device  are  known, 
functions  can  be  generated  in  such  a  way  that  these  limits  will  not  be  exceeded  during  the 
course  of  the  testing.  Also,  knowing  what  these  limits  are  will  enable  users  of  the  RWE  to 
anticipate  its  effectiveness  prior  to  a  test  to  allow  optimum  utilization  of  its  capabilities. 

In  order  to  calculate  the  maximum  angular  acceleration  of  the  RWE,  the  mass  moment 
of  inertia  of  each  rotating  component  must  be  determined.  The  first  component  examined 
is  one  of  the  rotating  louvers  which  is  shown  in  Figure  A-l. 
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DIAMOND  SECTION  J  CVUNORCAL  SECTION  -* 


NOTES: 

1.  NOT  TO  SCALE 

2.  ALL  DIMENSIONS  IN  INCHES 

3.  CL  -  2.830 

4.  T  -  0.3125 

5.  MATERIAL:  STEEL 

Figure  A-l:  Diamond  Shaped  Louver 


The  equation  used  to  calculate  the  mass  moment  of  inertia  for  the  diamond  shaped 
section  of  the  louvers  is  shown  in  Equation  (A-l). 


h  =  +  f.2) 


(A-l) 


In  this  equation,  p  is  the  density  of  steel  and  is  7849  kg/m3.  The  other  values  can  be  obtained 
from  Figure  A-l.  Solving  this  equation  yields  a  mass  moment  of  inertia  of  1.239  x  10-4  kg—m 2 
for  the  diamond  shaped  section  of  one  louver. 


The  inertia  of  the  remaining  portions  of  the  louver  must  be  calculated  next.  These 
two  sections  are  both  cylindrical  and  can  therefore  be  lumped  together  for  this  calculation. 
Equation  (A-2)  illustrates  the  means  of  calculating  this  value. 

pirlr 4 


L  = 


(A-2) 
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In  this  equation  /  is  the  total  length  of  the  cylindrical  sections  and  r  is  the  radius  of  the 
cylindrical  section.  Substitution  of  these  values  into  Equation  ( A-2)  yields  2.525  x  10-7  kg  — 
m2.  When  added  to  the  value  obtained  from  Equation  (A-l)  the  total  mass  moment  of 
inertia  for  one  louver  can  be  determined  as  seen  in  Equation  (A-3). 

Iiouver  =  Id  +  h  =  1.242  x  10"4  kg  -  m2  (A-3) 

Connected  to  each  louver  is  an  aluminum  gear  (p  =  2730  kg/m3).  The  mass  moment  of 
inertia  of  these  gears  must  also  be  calculated.  A  diagram  of  the  gear  is  given  in  Figure  A-2. 


NOTES: 


1.  ALL  DIMENSIONS  IN  INCHES 

2.  MATERIAL:  ALUMINUM 

Figure  A-2:  2.5  Inch  Louver  Gear 


The  mass  moment  of  inertia  for  this  gear  is  calculated  using  Equation  (A-4). 


(A-4) 


The  total  inertia  for  the  gear  is  determined  by  summing  the  inertia  values  calculated 
for  regions  A  and  B.  The  result  for  one  gear  is  2.782  x  10~5  kg  —  m2. 

From  this,  the  total  mass  moment  of  inertia  for  one  louver/gear  arrangement  is  calcu¬ 
lated  through  use  of  Equation  (A-5). 

//  =  I g^r  +  hower  =  1-520  x  10~4  kg  -  m2  (A~5) 


In  order  to  calculate  the  maximum  louver  angular  acceleration,  other  characteristics  of 
the  device  are  required.  These  are  summarized  below. 
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•  Peak  torque  of  the  motor  =  Tp  =  6.36  N  —  m 

•  Rotor  and  drive  shaft  inertia  =  Id  =  350.0  x  10-6  kg  —  m2 

•  Radius  of  louver  gear  =  rt  =  0.03175  m 

•  Radius  of  drive  gear  —  rd  =  0.01270  m 

•  Angular  acceleration  of  louvers  =  9t 

•  Angular  acceleration  of  drive  shaft  =  9d 

The  sum  of  torques  on  a  body  is  equal  to  the  mass  moment  of  inertia  of  the  body  about 
the  axis  of  rotation  multiplied  by  the  angular  acceleration  of  the  body.  By  performing  this 
analysis  on  the  gear  train  of  the  RWE,  the  maximum  angular  acceleration  can  be  obtained. 
Equation  (A-6)  is  the  summing  of  torques  on  the  drive  shaft. 

Tp  —  Ti  =  ld9d  (A-6) 


Since  the  four  louver /gear  assemblies  are  all  identical,  the  total  torque  that  they  provide 
to  the  drive  shaft  is  simply  four  times  that  of  one  louver  assembly.  This  torque  is  defined  in 
Equation  (A- 7). 

T,  =  (A— 7) 

n 

The  relationship  between  the  motion  of  the  drive  shaft  and  the  motion  of  a  louver  is  repre¬ 
sented  by  rd9d  =  r/0/.  By  substituting  this  equation  and  Equation  (A-7)  into  Equation  (A-6) 
the  maximum  louver  angular  acceleration  reduces  to  Equation  (A-8). 


[( Vi/rj )  -f  (4 hrd/ri)} 


(A-8) 


When  Equation  (A-8)  is  solved,  the  maximum  louver  angular  acceleration  turns  out  to  be 
5687.6  rad/s2. 

To  be  sure  that  the  capabilities  of  the  motor  would  not  be  exceeded,  the  maximum 
angular  acceleration  was  reduced  to  5000  rad/s2  from  5687.6  rad/s2. 

The  rated  speed  of  the  motor  is  listed  as  1100  rev/ min  which  is  equivalent  to  115.2  rad/s. 
By  using  the  relation  rd9d  —  rfii,  the  maximum  louver  angular  velocity  turns  out  to  be 
46.08  rad/s.  Again,  to  remain  safely  within  the  limitations  of  the  motor,  this  value  was 
reduced  to  43  rad/s.  By  using  Equation  (A-9),  the  number  of  motor  angular  units  (MAU) 
displaced  at  this  velocity  in  a  2  ms  time  interval  can  be  found. 


43  rad 


)  o  ( 


25000  MAU 
360° 


)  (0.002  s)  =  342  MAU  (A-9) 


With  these  values  known,  the  functions  can  be  developed  in  such  a  way  that  the  operational 
limits  of  the  motor  will  not  be  exceeded. 
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APPENDIX  B:  Pressure  Histories  for  all  Measurement  Stations 
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Toot  Static*  UoHm  ■  7J1I  motan  from  Mophrogm 


Toot  Station  LomNm  7.900  imtiri  from  Waphwgw 


Toot  Staff**  ImbHou  ■  8.281  motoro  from  fllophroQm 


Figure  B-l:  Pressure  Histories  for  9018  cm3  Driver;  Pressure  Level  1 
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Figure  B-2:  Pressure  Histories  for  9018  cm3  Driver;  Pressure  Level  2 
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Toot  Station  Location  ■  7.31 1  motor »  from  tXophm^m 


Too*  Station  Location  «  7.900  motor*  from  CHoptirogm 


»  i  i - t - i - r - 1 - 1  i - r - »  i  i 

-.01  OjOO  OlOI  042  003  004  006  006  007  006  OOf  OlO 

TJmo  After  Shack  Arrival  ($) 


Figure  B-3:  Pressure  Histories  for  9018  cm3  Driver;  Pressure  Level  3 
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Tmt  Stotfe*  Un«M  »  7J1»  m *an  from  Ot ophrotm 


Taaf  SlaHofl  Loaorton  ■  7.900  matara  from  Mophrogm 


Toaf  Station  laaoNon  »  ft. 221  mataro  tram  Moftftmgm 


Figure  B-4:  Pressure  Histories  for  9018  cm3  Driver;  Pressure  Level  4 
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Tima  After  Sham  Arrival  (a) 


Figure  B-5:  Pressure  Histories  for  4988  cm3  Driver;  Pressure  Level  1 
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root  Simla.  Loopkon  a  Mil  mNn  (run  Waptiragin 


Taat  3MIm  Losatton  a  7.900  m otoro  (run  OtQQtwaga. 


TVra  Attar  Shock  Arrival  («) 


Toot  Station  Location  a  Utl  motcro  from  Woohrogm 


Figure  B-6-  Pressure  Histories  for  4988  cm3  Driver;  Pressure  Level  2 
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-.01  OlOO  0.01  0.02  0X3  OO*  0X0  0X0 

Tim#  Attar  Shock  Arrival  (•) 


Tact  Station  Looaflon  »  7.900  mater*  tram  Diaphragm 


-.01  0.00  0X1  0X2  0X3  0X4  0X0  0X0 

Tima  Attar  Shock  Antral  (•) 


Figure  B-7:  Pressure  Histories  for  4988  cm3  Driver;  Pressure  Level  3 
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-.01  000  061  002  063  06*  009  006  007 

Tbiw  Af Hr  Stock  ArrM  («) 


Figure  B-8:  Pressure  Histories  for  4988  cm3  Driver;  Pressure  Level  4 
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SMI.  OMtawn  (hf.)  <•*•> 


Tee*  StotteA  lo«oK«n  ■  7411  meter*  from  Oteptetigm 


Figure  B-9:  Pressure  Histories  for  2840  cm3 


Driver; 


Pressure  Level  1 
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„un  <»«)  Stoll.  Owpmun  (kfo)  »««.  »»rp«..ur.  (kl>.) 


Toot  Slotted  loodtten  ■  7J1I  moNro  from  Diaphragm 


Toot  Station  Looatten  •  7.900  motero  from  Olophrogm 


Toot  Slot* on  Loootton  ■  A. 281  moloro  from  Dtephragm 


Figure  B-10:  Pressure  Histories  for  2840  cm3  Driver;  Pressure  Level  2 
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Toot  Station  Leoatten  ■  7.31*  motor*  Inn  Ptoptwagm 


Toot  Station  Loaatton  ■  8.281  motor*  from  Maptvagm 


Tim*  Aft*r  Shook  Arrt v*l  <•) 


Figure  B-ll:  Pressure  Histories  for  2840  cm3  Driver;  Pressure  Level  3 
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Toot  Hottoo  liogWm  *  7.51 1  motors  from  Diaphragm 


Toot  Station  LooaHon  •  7.900  motors  from  Ofophrogm 


Toot  Station  Location  ■  1211  motors  from  Otophrogn 


Figure  B- 12:  Pressure  Histories  for  2840  cm3  Driver;  Pressure  Level  4 
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ATTN:  Aerospace  Library 

P.O.  Box  8757,  BWI  Airport 

P.O.  Box  3707 

Baltimore,  MD  21240 

Seattle,  WA  98124 

1 

Director 

1 

California  Research  &  Technology,  Inc 

NASA-Langley  Research  Center 

ATTN:  M.  Rosenblatt 

ATTN:  Tech  Lib 

20943  Devonshire  Street 

Hampton,  VA  23665 

Chatsworth,  CA  91311 
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1  Carpenter  Research  Corporation 

ATTN:  H.  Jerry  Carpenter 
27520  Hawthorne  Blvd.,  Suite  263 
P.O.  Box  2490 

Rolling  Hills  Estates,  CA  90274 

1  Dynamics  Technology,  Inc. 

ATTN:  D.  T.  Hove 
Suite  300 

21311  Hawthorne  Blvd. 

Torrance,  CA  90503 

1  EATON  Corporation 
Defense  Valve  &  Actuator  Div. 
ATTN:  Mr.  J.  Wada 

2338  Alaska  Ave. 

El  Segundo,  CA  90245-4896 

2  FMC  Corporation 
Advanced  Systems  Center 
ATTN:  Mr.  J.  Drotleff 

Ms.  C.  Krebs,  MDP95 
2890  De  La  Cruz  Boulevard 
Box  58123 

Santa  Clara,  CA  95052 

1  Goodyear  Aerospace  Corporation 

ATTN:  R.  M.  Brown,  Bldg  1 
Shelter  Engineering 
Litchfield  Park,  AZ  85340 

4  Kaman  AviDyne 

ATTN:  Dr.  R.  Ruetenik  (2  cys) 
Mr.  S.  Criscione 
Mr.  R.  Milligan 
83  Second  Avenue 
Northwest  Industrial  Park 
Burlington,  MA  01830 


No  of 

Copies  Organization 

3  Kaman  Sciences  Corporation 
ATTN:  Library 

P.  A.  Ellis 
F.  H.  Shelton 

1500  Garden  of  the  Gods  Road 
Colorado  Springs,  CO  80907 

1  Kaman  Sciences  Corporation 
ATTN:  Mr.  F.  W.  Balicki 

6400  Uptown  Boulevard  N.E.  Suite  300 
Albuquerque,  NM  87110 

2  Kaman-TEMPO 
ATTN:  DASIAC 

Don  Sachs 
P.O.  Drawer  QQ 
816  State  Street 
Santa  Barbara,  CA  93102 

1  Ktech  Corporation 

ATTN:  Dr.  E.  Gaffney 

901  Pennsylvania  Avenue,  N.E. 

Albuquerque,  NM  87111 

1  Lockheed  Missiles  &  Space  Co. 

ATTN:  J.  J.  Murphy, 

Dept.  81-11,  Bldg.  154 
P.O.  Box  504 
Sunnyvale,  CA  94086 

2  McDonnell  Douglas  Astronautics 

Corporation 

ATTN:  Robert  W.  Halprin 
K.A.  Heinly 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 

1  Orlando  Technology,  Inc. 

ATTN:  Mr.  D.  Matuska 
P.O.  Box  855 
Shalimar,  FL  32579 
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Physics  International  Corporation 
2700  Merced  Street 
San  Leandro,  CA  94577 

R&D  Associates 
ATTN:  Technical  Library 
Dr.  Allan  Kuhl 
P.O.  Box  9695 
Marina  Del  Rey,  CA  90291 

R&D  Associates 
ATTN:  G-P.  Ganong 
P.O.  Box  9377 
Albuquerque,  NM  87119 

Science  Applications,  Inc. 

ATTN:  W.  Layson 

John  Cockayne 
P.O.  Box  1303 
1710  Goodridge  Drive 
McLean,  VA  22102 

Science  Applications  International 
Corp. 

ATTN:  Mr.  J.  Guest 
2109  Air  Park  Rd,  SE 
Albuquerque,  NM  87106 

Sparta,  Inc. 

Los  Angeles  Operations 
ATTN:  I.  B.  Osofsky 
3440  Carson  Street 
Torrance,  CA  90503 

Sunburst  Recovery,  Inc. 

ATTN:  Dr.  C.  Young 
P.O.  Box  2129 

Steamboat  Springs,  CO  80477 

Sverdrup  Technology,  Inc. 

ATTN:  R.  F.  Starr 
P.O.  Box  884 
Tullahoma,  TN  37388 


3  SRI  International 

ATTN:  Dr.  G.  R.  Abrahamson 
Dr.  J.  Gran 
Dr.  B.  Holmes 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

2  S-CUBED 

A  Division  of  Maxwell 
Laboratories,  Inc. 

ATTN:  C.  E.  Needham 

Dr.  Lynn  Kennedy 
2501  Yale  Blvd.,  SE 
Albuquerque,  NM  87106 

3  S-CUBED 

A  Division  of  Maxwell 
Laboratories,  Inc. 

ATTN:  Technical  Library 
R.  Duff 
K.  Pyatt 
PO  Box  1620 
La  Jolla,  CA  92037-1620 

1  Texas  Engineering  Experiment 

Station 

ATTN:  Dr.  D.  Anderson 

301  Engineering  Research  Center 

College  Station,  TX  77843 

1  Thermal  Science,  Inc. 

ATTN:  R.  Feldman 
2200  Cassens  Dr. 

St.  Louis,  MO  63026 

1  TRW 

Ballistic  Missile  Division 
ATTN:  H.  Korman, 

Mail  Station  526/614 
P.O.  Box  1310 
San  Bernadino,  CA  92402 
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1  Battelle  Memorial  Institute 

ATTN:  Technical  Library 
505  King  Avenue 
Columbus,  OH  43201 

1  California  Institute  of  Technology 
ATTN:  T.  J.  Ahrens 

1201  E.  California  Blvd. 

Pasadena,  CA  91109 

2  Denver  Research  Institute 
ATTN:  Mr.  J.  Wisotski 

Technical  Library 
P.O.  Box  10758 
Denver,  CO  80210 

1  Massachusetts  Institute  of  Technology 

Aeroelastic  and  Structures  Research 
Laboratory 

ATTN:  Dr.  E.  A.  Witmer 
Cambridge,  MA  02139 

1  Massachusetts  Institute  of  Technology 
ATTN:  Technical  Library 
Cambridge,  MA  02139 

2  New  Mexico  Engineering  Research 

Institute  (CERF) 

University  of  New  Mexico 
ATTN:  Dr.  J.  Leigh 
Dr.  R.  Newell 
P.O.  Box  25 

Albuquerque,  NM  87131 

1  Northrop  University 

ATTN:  Dr.  F.  B.  Safford 
5800  W.  Arbor  Vitae  St. 

Los  Angeles,  CA  90045 
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4  Southwest  Research  Institute 

ATTN:  Dr.  W.  E.  Baker 
A.  B.  Wenzel 
Dr.  C.  Anderson 
S.  Mullin 

6220  Culebra  Road 
San  Antonio,  TX  78284 

1  Stanford  University 

ATTN:  Dr.  D.  Bershader 
Durand  Laboratory 
Stanford,  CA  94305 

3  University  of  Maryland 

Department  of  Mechanical  Engineering 
ATTN:  Dr.  W.  Fourney 
Dr.  R.  Dick 
Dr.  J.  Williams 
College  Park,  MD  20742 

Aberdeen  Proving  Ground 
1  Cdr,  USATECOM 

ATTN:  AMSTE-TE-F  (L.  Teletski) 

1  Cdr,  USATHMA 

ATTN:  AMXTH-TE 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes. 
Your  comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  BRL  Report  Number  BRL-TR-3149 _ Date  of  Report  September  1990 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 

for  which  the  report  will  be  used.) _ 


1 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ 


j  . 

5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 

saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 

,  ...  ‘i _ • _ 

—————  - 

6. Gehdhd  Comments.  Whal  do  $du  think  should  be  changed  to  improve  future  reports?  (Indicate 

changes  *)  organization,  technical  content,  format,  etc.) _ _____ _ 


■  v'  ■*■,?*■ 

;v!7gggBBI - 

— ****i?t¥T —————— 

<»v 

■  "3 

Name 


CURRENT  Organization 

ADDRESS  _ 

Address 


City,  Slate,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Correct 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 

OLD 

ADDRESS 

Organization 

Address 

City,  State,  Zip  Code 

(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


FOLD  HERE 


Department  of  the  Army 

Director 

U.S.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21  (Mr  -5066 
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BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  No  0001,  APG,  MD 


POSTAGE  WILL  BE  PAID  BY  ADDRESSEE 


Director 

U.S.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21005-9989 


FOLD  HERE 


